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FOREWORD 

This document is Volume II of a two-volume report presenting the latest version of the 
Plume Impingement Computer Program (PUMP). The theoretical bases of the program 
and recent changes are discussed in Volume I. This volume describes the user input 
necessary to run the program and outlines two sample problems. 

The NASA/MSFC Contracting Officer’s Technical Representative for this contract was 
Mr. Peter Sulyma of the Induced Environments Branch, ED33, Aerophysics Division of 
the Structures and Dynamics Laboratory. 
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Section 1 
INTRODUCTION 

The Plume Impingement Program (PLIMP) is a computer code used to predict impact 
pressures, forces, moments, heating rates and contamination on surfaces due to direct 
impingement of flowfields. Typically, it has been used to analyze the effects of rocket 
exhaust plumes on nearby structures from ground level to the vacuum of space. The 
program normally uses flowfields generated by the MOC [1], RAMP2 [2], SPF/2 [3], or 
SFPGEN [4] computer programs. It is capable of analyzing gaseous and gas/particle 
flows. A number of simple subshapes are available to model the surfaces of any 
structure. 

The original PLIMP program [5] has been modified many times over the last 20 
years with Refs. [6,7] documenting the major changes. The theoretical bases for the 
referenced major changes, and additional undocumented changes and enhancements 
since 1988 are summarized in Volume I of this report. This volume is the User’s Input 
Guide and should be substituted for all previous guides when running the latest version of 
the program. This version can operate on VAX and UNIX machines with NCAR graphics 
ability. 


1 



REMTECH 


RTR 214-01 


Section 2 
INPUT GUIDE 

This section describes the input necessary to use the 1991 version of PUMP. 
Subsection 2.1 explains the input for the impingement calculations, Subsection 2.2 shows 
the files necessary to be opened and Subsection 2.3 explains the input for the new 
graphics program. 

2.1 Input Instructions 

RECORD 1 : PROGRAM CONTROL INFORMATION 
Format: (715,13,211,715) 

| Column j Variable* [ R Description 

5 IGO Type of calculations desired "~~™ 

= 0 Perform impingement analysis only 
= 1 Arrange data generated by RAMP2 or MOC program into a 
format used by PLIMP and perform an impingement 
analysis. 

= 2 Arrange data only 

10 I RUN (4) Source of impingement data — Used to indicate if 

impingement calculations are to be made or if impingement 
data are to be read into the program from a prior execution. 

= 0 Impingement data was not previously calculated 
= N Fortran file number containing calculated impingement data 
from a previous PLIMP run. 

NOTE: If IRUN(4) > 0 only Record 1, NAMELIST INPT2 (see after Record 1 9), and 

Record 20 are used. 

~ 15 IRUN(6) Not used """ 

20 IRUN(7) Output type (used if IRUN(4) > 0, Column 10) 

= 0 Do not print data 
= 1 Print data only 
= 2 Plot data only 
= 3 Print and plot data 
25 LSOLID Not used 

30 LRAMP Input flowfield type 

= 0 Any input flowfield file except RAMP2, MOC and SFPGEN 
= 1 RAMP2, MOC or SFPGEN input flowfield file 
35 IAF(8) Type of flowfield chemistry 

= 0 Variable O/F, constant O/F, two-phase, ideal gas, variable 
enthalpy 
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Column 

Variable 

Description | 


= 1 

Variable total enthalpy (single phase) 


= 2 

Ideal gas with variable total enthalpy 

38 

ISPECI 

Debug flag for subroutines SPECIE and SPLOAD - Species 
mole fractions, O/F ratios, entropy level 


= 0 

Normal output 


= 1 

Extra output 

39 

IDSHAD 

Debug flag for subroutine SHADOW - Shading of subshapes 


= 0 

Normal output 


= 1 

Extra output 

40 

IMTRAN 

Debug flag for subroutine MTRAN - Transitional region 
calculations 


= 0 

Normal output 


= 1 

Extra output 

45 

IFRTAP 

Not used 

50 

ISORCG 

Gaseous source flowfield 


= 0 

No gaseous source flowfield 


= 1 

Gaseous source flowfield 

55 

ISORCP 

Particle source flowfield 


= 0 

No particle source flowfield 


= 1 

Particle source flowfield 

60 

METRIC 

Type of units 


= 0 

Input and output files (units 5 and 6) English units 


- 1 

Input and output files (units 5 and 6) SI units 

NOTE: 

Debug flags will output In English units only. 1 

65 

FORM 

Output format 


= 0 

Output file (unit 6) in 132 column format 


= 1 

Output file (unit 6) in 80 column format 

70 

NPLOT 

NCAR graphics data file for pre- or postprocessing 


= 0 

Do not print impingement data to unit 9 


= 1 

Print impingement data to unit 9 

75 

ICENT 

Not used if NPLOT = 0 


= 0 

Unit 9 data for postprocessor 


= 1 

Preprocessor data (unit 9) do not run impingement analysis 
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RECORD 2: CONTROL INFORMATION FOR THE ORDERING SECTION (used only 
when IGO > 0 Record 1, column 5) 


Format: 

(1615) 



Column 

1 Variable 

J_ Description 

5 

ISTART 

= N 

Beginning of flowfield ordering calculations 
Ordering the flowfield data will begin with this 
characteristic line number (MOC) or normal data surface 
(RAMP2) 

10 

ISIGN 

= -1 

Location of flowfield quadrant (MOC runs) 

The flowfield was generated in the fourth quadrant 



= 0 

The flowfield was generated in the first quadrant 

15 

KNUMBR 


Not used 

20 

IDEL 

= 0 

Deletion of flowfield points 

Will not delete any flowfield points before ordering 



= 1 

Will delete flowfield points (based on Record 3, Columns 
21-30) 

25 

IPRINT 

= 0 

Printout options 

No intermediate printout while ordering data 



= 1 

Intermediate data are to be printed as the flowfield is 
ordered by distance from the engine exit plane 

30 

ITERM 

- N 

End of flowfield ordering calculations 
Characteristic line or normal number where ordering of 
flowfield data is to be terminated. [If 0, it is set to 
1 ,000,000 . . . flowfield calculations (MOC, RAMP2) 
terminated properly] 

35 

ISEND 

= 1 

Boundary determination 

Plume boundary is to be curve-fitted for use in the 
interpolation scheme 



= 2 

Cutoff limits read as input data (Record 18) are to be 
used to see if a point is within a prescribed boundary 




(this option is usually used). 

40 

ISKIP 


Deletion of lines in boundary determination (used if 
ISEND = 1) 



= 0 

Every line will be examined for a free boundary point 



= N 

Every Nth line will be examined for a free boundary point 

45 

ISKPL 

= 0 

Deletion of characteristic lines in ordering process 
Read all characteristics 



= N 

Read every Nth + 1 characteristic 
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RECORD 3: REFERENCE COORDINATES FOR THE ORDERING CALCULATIONS 
(used only when IGO > 0 Record 1, column 5) 

Format: (6E10.6) 


Column 

Variable 

Description 

1-10 

RREF 

Radial coordinate to which each flowfield data point will 
be referenced 

11-20 

XREF 

Axial coordinate to which each flowfield data point will 
be referenced 

NOTE: 

XREF and RREF are used to accomplish any desired coordinate system 
translation. This provides a simple means to move the origin; e.g., from the 
nozzle throat to the nozzle exit.) [X=(X-XREF)/DIAM, R=(R-RREF)/DIAM] 

21-30 

DELETE 

One of two points with a distance between them less 
than DELETE will be deleted from the flowfield data 
(used if IDEL = 0, Record 2, column 20) 

31-40 

DIAM > 0.0 

Reference factor, units consistent with plume 

dimensions. Can be used to scale, etc., the local plume 

coordinates 

(use 1 .0 for no scaling) 
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RECORD 4: RUN IDENTIFICATION 
Format: (18A4) 


Column : Variable Description 

1-72 HEADER Title or heading 



Figure 1 : Elemental Strips 
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RECORD 5: INPUT OPTIONS NECESSARY TO CONTROL PROGRAM EXECUTION 
Format: (1615) 


Column Variable Description 


5 

IOPT(1) 

= 0 

Flowfield type 

Flowfield is two-dimensional 



= 1 

Flowfield is axisymmetric 

10 

IOPT(2) 

= 0 

Desired heating options 

Calculate continuum and free molecular heating 



= 1 

Calculate continuum heating only 



= 2 

Calculate free molecular heating only 



= 3 

None 

15 

IOPT(3) 

= 0 

Uniform flow assumptions (usually use 3) 
Composite vehicle in uniform flow 



= 1 

Subshapes in uniform flow 



= 2 

Elemental strips in uniform flow. (See Fig.1, page 6.) 



= 3 

Elemental areas in uniform flow 

19-20 

IOPT(4) 

= N 

Number of subshapes in composite body (maximum of 
50 subshapes) 

22 

IVAPN 

= N 

Number of species whose evaporation rate data as a 
function of surface temperature is to be input (5 max) [At 
present, data for N 2 H 4 , H 2 0, NH 3> C0 2 and HCI are 
included. Use this option only if you want to add species 
and ISCTAM * 2 (this record, column 24).] 

23 

IPROP 


Indicates motor propellant type [Used only when net 
deposition rate output (ISCTAM=2) is selected.] 



= 0 

ISCTAM < 2 



= 1 

Hydrazine - cold start 



= 2 

Hydrazine - pulsed mode 



= 3 

MMH/N 2 0 4 bipropellant 



= 4 

Solid propellant 



= 5 

Other 

24 

ISCTAM 


Impingement gaseous contaminate species calculation 
(contaminant species that are impinging on the surfaces) 



= 0 

Do not calculate and output species at surface 



= 1 

Calculate and print out species which exist at the surface 
and are available to deposit on surface 



= 2 

Calculate net deposition rates using model described in 
Section 3.3.4, Vol. 1. 
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Column Variable 

Description 

25 IOPT(5) 

Type output desired 

= 0 

Full output (data at each elemental area) 

= 1 

Force and torque summary only 

30 IOPT(6) 

Local flowfield chemistry 

= 0 

Consider property variation with entropy 

= 1 

Neglect entropy effects (use only one gas table) 

35 IOPT(7) 

Boundary layer chemical species 

= 0 

Use species on input flowfield file 

= N 

Number of species input (requires additional data from 


Record 10) (25 maximum) 

37 IOPT(8) 

Forces and moment scaling (usually use 0) 

= 0 

Do not scale forces and moment 

= 1 

Multiply resultant forces and moment by PC RAT (Record 


10, columns 61-70) 

38 =0 

Plume loading only 

= 1 

Subtract P Noo (Newtonian) from P impaci for points within 
plume 

= 2 

Calculate PN Ximpact lor points outside plume 

= 3 

Subtract P^ M over entire body surface 

39 =0 

Zero pressure used for points in plume but shaded from 

flow 

= 1 

External P^ used for points in plume boundary but 

shaded from the flow 

O 

11 

o 

Single case being run 

= 1 

Multiple cases with limited changes to input data (new 

cases input with NAMELIST CHNGE after Record 19). 

45 IOPT(9) 

Type of local impact pressure calculation (continuum 

regime only, Section 3.2.3, Vol. 1) 

= 1 

Newtonian assumption ( K cp = 2.0) 

= 2 

Modified Newtonian [K cp = (0.814 + ) C ps \ 

= 3 

sentropic expansion 

= 4 Oblique shock 

- 5 

Limited modified Newtonian ( K cp = C ps ) 

NOTE: Record 14 provides opportunity to select a particular calculation method for 

1 each subshape. 
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Column Variable 

• ■ - Description | 

46-50 

IOPT(10) 

= 0 

Debug Flag. 

No debug printout 



= N 

Number of desired subshapes in debug printout. (Look 
out! starts with subshape 1 and continues through N 




subshapes.) 

51-55 

IOPT(1 1) 


Plot control for SC 4020 plot program (Not NCAR 
graphics) 

NOTE: 

It may be necessary for user to update PUMP FORTRAN file to output 
parameters in the format necessary for his particular Dlot Droaram. 

51 


= 0 

Do not connect plotted data points 



= 1 

Connect plotted data points 

52-53 


= N 

Number of file unit on which impingement plot data will 
be written, usually 14 

54-55 

ITPLOT 

= N 

Number of parameters to be plotted 

56-60 

IOPT(12) 

= N 

Case number 

65 

IOPT(13) 

= 0 

Dimensions of input 

Impingement program dimensions in inches 



= 1 

Impingement program dimensions in feet (meters if 
METRIC = 1, Record 1, column 60) 

70 

IOPT(14) 

= 0 

Dimensions of ordered flowfield 
Flowfield dimensions in inches 



= 1 

Flowfield dimensions in feet 

75 

IOPT(15) 

= 0 

Local property printout 
No intermediate printout of properties 



= 1 

Print local properties 



= 8 

Intermediate printouts in heating routines 

80 

IOPT(16) 

= 0 

Body-on-body shading 

No body-on-body shading calculations to be performed 



= 1 

Body-on-body surface shading to be considered. (Input 
additional shading data after each subshape geometry, 




Record 14c) 


9 




f= 3 rvxi -T* e: g »— i 


RTR 214-01 


RECORD 6: REFERENCE COORDINATE SYSTEM POSITION VECTOR (Position 
vector defining the relative positions of the engine and reference coordinate systems.) 

This record contains values for an engine position vector, REC, a moment transfer 
vector, RECT, and three control flags, ITU, ITV, ITW. These flags control the use of 
the position vector and its transformation matrix, and the moment transfer vector. REC 
and RECT are measured in inches, feet, or meters consistent with IOPT(13) (Record 
5, column 65). 


yR 



Engine coordinates may be 
defined from reference sys- 
tem, or reference system 
may be defined by engine 
system. 


Sketch of Coordinate System Convention 

NOTE: Rotations are always from the reference coordinate system toward the engine 
coordinate system. All coordinate rotations discussed in this report obey the 
right-hand rule for positive rotation notation. See sketch below. 


Z R 


Z 


t 


yR,y' 





Right-Hand Convention for Positive Angles of Rotation 
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Format: ( 6E1 0.6,31 5) 


Column 

Variable 

Description 

1-10 

REC(x) 


11-20 

REC(y) 

Components of engine position vector 

21-30 

REC(z) 


31-40 

RECT(x) 

Components of moment transfer vector written in reference 

41-50 

RECT(v) 

51-60 

RECT(z) 

system (RECT = 0 if ITW = 0, column 75) 

65 

ITU 

Engine position vector control flag 


= 0 

The engine position vector, REC, is measured in the engine 
system and locates the reference with respect to the engine. 
(Vector points from the engine to the reference system.) 


= 1 

The position vector, REC, is measured in the reference co- 
ordinate system and locates the engine system with respect 
to the reference system. (Vector points from the reference to 
the engine system.) 

68-70 

ITV 

Engine reference system transformation matrix control flag 
(Input is right adjusted.) 


= 0 

The components of the i and j unit vectors of the transforma- 
tion matrix will be read from Record 7, columns 1-60. 


>0 

The three axes angles of rotation will be read from Record 
7. The program will compute the transformation matrix. User 
must specify order of rotation(s) with this flag, ITV, as a three- 
digit coded number using the following (rotate reference to 
engine system): 


- 1 

x-axis rotation, <f> deg 


= 2 

y-axis rotation, 9 deg 


= 3 

z-axis rotation, if> deg 

75 

ITW 

Moment reference point control. 


NOTE: PUMP will normally calculate moments in the reference system* However, 
the reference system origin is typically chosen for convenient subshape 
description. It may not be the point about which the surface would rotate 
if a moment were applied. This option allows the user to define that point 
(e.g., the surface’s center of gravity). PUMP will now calculate all moments 
in this coordinate system. : • ' ■ ' 

= 0 The resultant moments are referenced to reference coordinate 
system 

= 1 The resultant moments are to be referenced to the point 
defined by RECT 
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RECORD 7: REFERENCE COORDINATE SYSTEM DIRECTION COSINES OR AXES 
ROTATION ANGLES 

This record contains values for the direction cosines of the reference system unit 
vectors, 1C and JC. The unit vectors are measured in the engine, E, coordinate system. 
If the control variable, ITV (Record 6, columns 68—70), is greater than zero, the axis 
rotation angles, <p, 9, ip, are input and the cosines computed by the program. 

Format: (6E10.6, 15) 


Column 1 

Variable 

Description 

1-10 

IC(x) 

Direction cosines of x R in E system 
than zero, then <p, 9, ip are input. 
Direction cosine from xe or <p (deg) 

if ITV=0. If ITV greater 
(Angles read in this order 

11-20 

IC(y) 

Direction cosine from y E or 9 (deg) 

regardless of the order of 

21-30 

IC(z) 

Direction cosine from z E or xp (deg) 

rotation(s), ITV) 

31-40 

UNITIC 

x, y, z direction cosines of y R in E system if ITV = 0. Not used 
if ITV greater than zero. 

Direction cosine from x E 

41-50 

UNITJC 

Direction cosine from y E 


51-60 

UNITKC 

Direction cosine from z E 


61-65 

IFL 

Must match ITV, Record 6, columns 68-70 
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RECORD 8: PLOT IDENTIFICATION PARAMETERS 
Format: (1615) 


If IOPT(11) > 0 (Record 5, columns 51-55) the following impingement parameters 
can be plotted by the SC4020 or user plotting routine inserting the correct number (right 
justified). (Again, the user will need to adapt his program to PLIMP format or vice versa.) 



1-5 IPLSYM(I) =1 Turbulent convective heating rate 


: = 2 Laminar convective heating rate 

= 3 Transitional flow convective heating rate 
66-70 IPLSYM(14) =4 Free molecular flow convective heating rate 

= 5 Local Reynolds number based on momentum thickness 
= 6 Impact pressure 

= 7 Impact pressure/chamber pressure ratio 
» 8 Impact force 
= 9 Incident mass flow rate 
= 10 Free-stream Mach number 
= 1 1 Local Mach number 
= 1 2 Local plume impingement angle 
= 1 3 Local temperature 
= 14 O/F ratio 


75 N Point deletion criterion. Every Nth output point will be 

plotted. 

80 N Line deletion criterion. Every Nth output line will be 

plotted. 


NOTE: The numbers input in columns 1-5 ... 66-70 do not have to be in : ^ 

numerically ascending order; however, these columns, in blocks of five, 
should be filled from left to right. Blanks should be used through 70, i.e, 
Column = (1 5)(6 1 0)(1 1 1 5)(1 6 20)... (61 65)(66 70)(71 75)(76 80) 

Number = 6 1 8 1 i 
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RECORD 9: SPECIES CONCENTRATION CARDS [used when IOPT(7) > 0 (Record 
5, column 35)]. 

This record specifies the mole fractions of the flowfield gas composition. It is input if 
the flowfield was generated by the MOC program using thermodynamic data input by file 
(ICON(1) = 1) or the RAMP2 program using the ideal gas and variable total enthalpy 
option (ICON(1) = 32). Input only the 25 most dominant species. Do not input if heating 
rate or contamination option is not chosen; i.e. p IOPT(2) = 3 or ISCTAM = 0 (Record 5, 
columns 10 and 24, respectively) or using SFPGEN or SPF/2 flowfields. 

Format: 3(2A4,2X,E10.6) 


Column 

1-8 


11-20 

21-28 

31-40 

41-48 

51-60 


Variable 

XPNAME(I) 


. ^ ; ■ iit i : ■: £ V: ; : : : :■ ; : ^ Description _ : 1: 0 , : V 

Chemical species name; i.e., H2, C02, etc. (See Table 1, 
page 15.) The species names are left justified. 


XP(1) 

XPNAME(2) 

XP(2) 

XPNAME(3) 

XPNAME(4) 


Mole fraction associated with chemical species 1 
Species name 2 

Mole fraction associated with chemical species 2 
Species name 3 

Mole fraction associated with chemical species 3 


Use additional records as necessary. 
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Table 1: Chemical Systems in PUMP 


*AI 

*cf 4 

*h 2 s 

Na 2 

AICI 

CH 

He 

*Ne 

AICIH 3 

ICH 4 

1 

0 

AIF 

CN 

*l 2 

OH 

AIF 3 

‘CO 

‘Li 

o 2 

AIN 

‘COS 

*LiCI 

P 

AIO 

* _ _ * * 
co 2 

*LiF 

PCI 3 

*Ar 

CP 

LiO 

PF 

B 

cs 

*Li 2 

ph 3 

BCI 

C 2 

Li 3 0 

PN 

bci 2 

c 2 h 2 

Mg 

PO 

*bci 3 

*c 2 h 4 

MgCI 

PS 

BF 

*c 2 n 2 

MgCI 2 

P 2 

bf 2 

Cl 

MgF 

P 4 

*bf 3 

CICN 

MgF2 

*s 

BO 

*CIF 

N 

*sf 6 

b 2 

*CIF 3 

*n 2 

so 

b 2 o 3 

‘CIO 

nf 3 

*so 2 

Be 

*CI 2 

NH 

s 2 

BeCI 

*F 

*nh 3 “ 

*Si 

BeCI 2 

FCN 

‘NO 

SiCI 

BeF 

*f 2 

‘NOCI 

SiCI 4 

BeF 2 

H 

‘N 2 0 

SiF 

*C 

*hcn“ 

‘Na 

SiF 4 

CCI 

*HCI“ 

*NaCI 

SiH 4 

CF 

HF 

‘NaF 

SiO 

cf 2 

*h 2 

NaO 

*S0 2 

cf 3 

*h 2 o" 

‘NaOH 

SiS 

Sl 2 

*UA (unreacting ammonia) 


*Vapor Pressure Data Available 
‘‘Evaporation Rate Data Available 
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RECORD 10: ENGINE REFERENCE AND NON-DIMENSIONALIZING PARAMETERS 
Format: (8E10.6) 


Column 

Variable 

Description 

1-10 

XTE 

Transfer distance from engine gimbal to exit plane (inches, 
feet, or meters consistent with IOPT(13), Record 5, column 
65). This option used primarily when the orientation of the 
engine coordinate system with respect to the reference 
system is about the engine gimbal point. 

11-20 

XMIN 

Minimum distance from engine reference system at which 
impingement calculations are to start. (This eliminates 
useless calculations.) 

21-30 


Not used 

31-40 

REFERN 

Engine scale factor if plume coordinates are 
non-dimensional, or to be scaled (units consistent with 
IOPT(14), Record 5, column 70). 

41-50 

RICH 

Diameter of nozzle at exit plane for use in 
non-dimensionalizing program output (units consistent with 
IOPT(13) Record 5, column 65). (Usually 1.0) 

51-60 

PC 

Engine chamber pressure (psi or pascals), 
non-dimensionalizes pressure output. (Input 1 .0 for no 
non-dimensionalization.) 

61-70 

PC RAT 

Ratio of desired plume chamber pressure for which heating 
rates, forces and moments are desired to the chamber 
pressure at which the plume flowfield was generated. 
(Usually 1 .0) 

71-80 

DELPHI 

Constant angular increment off stagnation line (deg) at which 
impingement calculations are desired. 
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RECORD 11: CONDITIONS EXTERNAL TO FLOWFIELD 


Format: 

(6E10.6) 


Column 

Variable 

} Description | 

1-10 

PAMB 

Static pressure, psfa, pascals 

11-20 

TAMB 

Static temperature, °R or °K 

21-30 

EMAMB 

Mach number 

31-40 

GAMAMB 

Ratio of specific heats, 7 

41-50 

RAMB 

Gas constant, f^/sec^-R or rr^/sec^-K 

51-60 

CPINF 

Newtonian pressure coefficient, C pa 
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RECORD 12 : HEATING RATE REFERENCE PARAMETERS 

This record must be input. [The record should be blank if no heating rates are to be 
calculated; i.e., IOPT( 2 ) = 3 (Record 5, column 10 )] 

Format: (7E10.6,1X,A4) 


Column 1 Variable I Description } 

I - 10 CPS Specific heat of solid particles BTU/lbm or kJ/kg(use 

0.33 for AI 2 O 3 ) 

II- 20 TAME Trajectory time if surface temperature calculations are to 

be made. 

21-30 XSTAG Axial distance from subshape origin to stagnation point. 

Units consistent with IOPT(13), (Record 5 , column 65) 
31-40 TD Temperature of the gas in the boundary layer at the 

stagnation point at onset of dissociation, °R or °K 
41-50 RB The reference radius (r*) at the stagnation point for use 

in velocity gradient calculation. (Units for r*. consistent 
with columns 51-60 this record.) 

51-60 FUDGE = 1 r* input in feet or meters 


= 12 r* input in inches 

(any numerical value input other than 1 or 12 defaults to 1) 
61-70 ACOMCF Accommodation coefficient for solid particle 

impingement ( 0.0 - 1 . 0 ) (usually 0 . 5 ) 

72-75 SHAP Nose shape of surface 

= 3DHM, Hemisphere 
= 2DHM, Cylinder 
= 3DDS, Circular disk 
= 2DFP, Rectangular flat plate 
= 3DCN, Cone 
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RECORD 1 3a: SUBSHAPE SURFACE TEMPERATURE [Input only if heating or contam- 
ination desired (IOPT(2) < 3 or ISCTAM > 0)] (Record 5, columns 10 and 24, respectively) 
This record contains the surface temperature of each of the N subshapes to be used 
in the heating rate calculation. N = IOPT(4) (Record 5, columns 19-20). 

Format: (8E10.4) 


I- 10 TWAL(1 ) Wall temperature (°R or °K) of subshape 1 

II- 20 TWAL(2) Wall temperature (°R or °K) of subshape 2 

21-30 TWAL(3) Wall temperature (°R or °K) of subshape 3 

• • 

• * 

71-80 TWAL(8) Wall temperature (°R or °K) of subshape 8 

Use additional records as necessary 


Variable 


Description 


Column 
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RECORD 13b: EVAPORATION RATE CONTROL (Input only if ISCTAM = 2 and 
IVAPN > 0, Record 5, column 22) 

Format: (2A4.I2) 


Column 

1-8 


15 


Variable 


AME 

IN 


: : |§|i Description ' : ■ v ; : ■;?•■■■" v”:# 

Species name (consistent with Table 1 and left adjusted). 

Number of evaporation rates vs. temperature points to be read 
in on Record 13c. 


RECORD 13c: SURFACE EVAPORATION RATES AS FUNCTION OF TEMPERATURE 
[Input data from lowest to highest temperature. Input only if IVAPN = 0 and ISCTAM = 
2 (Record 5, columns 22 and 24, respectively).] 

Format: (8E10.6) 


Column 

Variable 

Description 

EVPRAT(I) = T(1) + EVP(1) 

1-10 

T(1) 

Temperature - °R or °K — Lowest temperature point 

11-20 

EVP(1) 

Evaporation rate - lbm/ft 2 -sec or kg/m 2 -sec — Matches lowest 
temperature 

21-30 

T(2) 

Temperature - °R or °K — Next lowest temperature point 

31-40 

EVP(2) 

Evaporation rate - lbm/ft 2 -sec or kg/m 2 -sec — Matches next 
lowest temperature 


Continue until IN data sets have been input using additional records as 
necessary. 


NOTE: Records 13b and 13c are required only if the user has selected the deposition 
model (ISCTAM = 2, Record 5, column 2) and has indicated that there are ad- 
ditional surface evaporation rate data for species that are not already included 
in the internal program data. Input as many sets of Records 13b, c as IVAPN 
(maximum of 5). 
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RECORD 14a: REFERENCE DIMENSIONS 

This record contains N subshape reference dimensions to be used in determining 
the local flow regime. N = IOPT(4) p Record 5, columns 19-20. 

Format: (6E10.6) 

I Column | Variable 1 Description 

I- 10 REFDLT(1,1) Subshape reference length (inches, feet, or meters 

consistent with IOPT(13), Record 5, column 65). 

II- 20 REFDLT(2,1) Subshape 1 reference diameter (inches, feet, or meters) 
21-30 REFDLT(1,2) Subshape 2 reference length 
31-40 REFDLT(2.2) Subshape 2 reference diameter 
41-50 REFDLT(1,3) Subshape 3 reference length 
51-60 REFDLT(2,3) Subshape 3 reference diameter 

Use additional records as necessary. 
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RECORD 14b: SUBSHAPE GEOMETRY CARD 

These records identify the type of surface (conic, rectangular plate, circular plate, 
etc.) and the coefficients of the surface equation. There will be as many of these 
records as there are subshapes (maximum of 10). Records 14b, 14c, 14d are read as 
a set for each subshape. 

Format: (II ,I9,7E1 0.6) 


Column 


1 


Variable | Description 


IDSURF Surface identification 
= 1 Conic 

= 2 Rectangular plate 
= 3 Circular plate 
= 4 Right triangular plate 
= 5 Not used 

= 6 Body of revolution described by a polynomial 
= 7 Trapezoid 

NPHI = N Surface integration control 

NX = N Surface integration control. These are applied as follows: 


Subshape 

1. Conic 

2. R. Plate 

3. C. Plate 

4. T. Plate 

5. Not used 

6. Polynomial 

7. Trapezoid 


NPHI 

Angular (60 max) 

Z (unlimited) 
Angular (unlimited) 
Z (same as Y) 

Angular (60 max) 

Z (unlimited 


NX 

X (unlimited) 

Y (80 maximum) 
Radial (60 max) 

Y (60 max) 

X (unlimited) 

Y (60 max 


NOTE; The NPHI increment applies to the rectangular portion of the trapezoidal 
plate. Each triangular portion of the plate will have the same number of Z 
increments as there are Y increments. 


8 * N The type of impingement pressure calculation for the 

subshape if different from option specified by IOPT(9) 
(Record 5, column 45). If left blank, the option specified by 
IOPT(9) will be used. 

9 =0 Calculate heating rates on all surface strips (Fig.1, page 6) 
= 1 Heating rate calculations or data output on specified 

subshape strip numbers, if IOPT(5) = 1 (Record 5, column 
25), strip numbers specified by Record 14d. 

= 2 No heating rate calculations on subshape 
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Column Variable 


Description 


10 IHT 


= 3 Surface temperature calculations on selected strips specified 
by Record 1 4d. 

Continuum regime heat transfer rate method 
= 1 Eckert Reference Enthalpy method for unyawed cones 
= 2 Eckert Reference Enthalpy method for flat plates 
= 3 Yawed cone heating rates using swept infinite cylinder theory 
= 4 Stagnation point heating rates using Marvin and Deiwert 
Method (requires input data on Record 12) 

= 5 Cylinder heating rates using yawed infinite cylinder theory 


11-20 A 
21-30 B 
31-40 C 
41-50 D 
51-60 E 
61-70 XN 
71-80 XB 


Coefficients of subshape equations, subshape coordinate 
system. Refer to sketches below for definitions of 
coefficients; units of coefficients are consistent with IOPT(13), 
Record 5, column 65. Subshape coordinate system 
orientations must be located as shown. 


NOTE: For flat plates (rectangular, triangular, circular, trapezoidal), the X axis must 
be oriented in a general direction toward the engine. (Impingement assumed 
on side with +X axis.) 
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Conic , IDSURF = 1 



Rectangular Plate . IDSURF = 2 



Surface Equation: 

R = A (yB + CX + WO + E ), for XN <X< XB 
(E must = 0 for proper shading calculation.) 


Coordinate system origin is located at the geometric 
center. The X axis is normal to the plane of the plate. 
The Y axis should be oriented in the general direction 
of flow if heating rate calculations are made. 


A — 2Y mox 
B = 2Z moI 

(C, D, E, XN, XB unused.) 


Circular Plate. IDSURF = 3 



IDSURF - 4 



Coordinate system origin is located at the geometric 
center. The X axis is normal to the plane of the 
circular plate. Orientation of the Y, Z axes directions 
are established by the user. 

A = outer radius of circular plate, R<, 

B = inner radius of circular plate, Ri 
(C, D, E, XN, XB unused.) 


Coordinate system origin is located at the vertex of 
the right angle. The X axis is normal to plane of 
the plate. Orientation of the Y,Z axis should be in 
general direction of flow. 


A — Y max 

B = Zmax 

(C, D, E, XN, XB unused.) 
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Body of Revolution , IDSURF = 6 



Surface Equation 

R = A + BX + CX 2 + DX 3 + EX \ XN <X <XB 


Trapezoidal Plate . IDSURF = 7 


b - — C 



The coordinate system origin is at the base in the 
middle of the rectangle. The program will work from 
the +Z toward — Z direction. (The trapezoid must be 
symmetric about the Y-zxis.) 


A — V 

n ~ 1 max 

B = Base length 

C = Top length 

D = 1 Flow is from +X direction (impinges on 

front side) 

= —1 Row is from —X direction (impinges on 
back side) 

E = 0 Complete trapezoid 
= 1 +Z side of trapezoid 
= -1 -Z side of trapezoid 
= 2 Triangle on +Z side only 
= -2 Triangle on -Z side only 
= -3 Rectangular part only 
(XN and XB unused.) 


NOTE: The trapezoid routine can analyze a triangular plate and a rectangular plate; 
however, a substantial amount of execution time may be conserved if these 
simple subshapes are analyzed with their respective routines. The reason for 
this is due to the more complex computer logic involved in the trapezoid routine; 
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RECORD 14c: BODY SHADING DATA (input only if shading option is flagged, 
IOPT(16) = 1, Record 5, column 80) 

Format: (1012) 


Column 

Variable 

Description 

1-2 

3-4 

5-6 

19-20 

ICHEK(I) 

ICHEK(2) 

ICHEK(3) 

ICHEK(IO) 

Relative number(s) of the body(s) according to the order in 
which they are input which are to be checked for its shading 
on the body currently being input. (Input body numbers right 
justified.) 

NOTE: 

A maximum of 10 surface numbers may be Input. Scan for surfaces to check 
is terminated by reading a zero if less than 10 surfaces are Input. 

RECORD 14d: SURFACE TEMPERATURE IDENTIFICATION (input only if flagged on 
Record 14b, column 9 = 1 or 3) This option allows a simple thermal analysis of five 
points on a surface so that the surface temperature can vary with time. 

Format: (5I6.5E1 0.6) 

Column 

Variable 

lilili Description iPilM 


1-6 

ISELT(I) 

First point 

1-2 

N 

Strip number for first point (Fig.1, page 6) 

3-4 

5-6 

N 

= 1 
= 2 
= 3 

Node number on strip for first point (node numbers start at 
Ymin) 

Surface material 
Titanium 
Carbon/Carbon 
Aluminum 

7-12 

ISELT(2) 

Second point 

25-30 

ISELT(5) 

Repeat until all points are identified (5 maximum) 

31-40 

THK(1) 

Material thickness [units consistent with IOPT(13) Record 5, 
column 65] for first point 

41-50 

THK(2) 

Material thickness for second point 

71-80 

THK(5) 

Repeat until material thicknesses for all points are read in. 
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RECORD 15: SUBSHAPE POSITION VECTOR 

These records contain the position vectors, RCI, from the reference system to the 
subshape origins. RCI is measured in the reference system and “points" to the subshape 
origin. There will be as many of these records as there are subshapes which are read 
successively. The units are consistent with IOPT(13), Record 5, column 65. These 
records are read in the same order as the subshape geometry data (Record 14). 
Format: (3E1 0.6,15) 


Column 

Variable 

. ■ Description : : ■ ; 

1-10 

RCI(x) 

x component of position vector 

11-20 

RCI(y) 

y component of position vector 

21-30 

RCI(z) 

z component of position vector 

33-35 

ITA 

= 0 
> 0 

Reference-Subshape system status flag 

Matrix of direction cosines will be read 

Angles of axes rotations are to be read. Program will 

compute transformation matrix. User must specify order 

of rotation(s) using the following code: 

= 1 Rotation about x axis, <f> 

= 2 Rotation about y axis, 9 deg 
= 3 Rotation about z axis, ^ deg 

NOTE: 

The rotation is at 

ways from the reference to the subshape coordinate system. 
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RECORD 16: SUBSHAPE DIRECTION COSINES OR AXES ROTATION ANGLES 
These records contain either the subshape direction cosines of the unit vectors, II 
and Jl, or the axes angles of rotation, depending on the value of ITA (Record 15, columns 
33-35). There will be one of these records for each subshape. These parameters are 
measured in the composite system reference frame. If the control variable ITA is greater 
than zero, the axis rotation angles, <p, 9, ip, are input and the cosines computed by the 
program. These records are read in the same order as the geometry data (Record 14) 
immediately after Record 15. 

Format: (6E1 0.6,15) 


Column 


Variable 


: 1 1 j 1 1 III : 1 | II Description- s : i:- ; 

Direction cosines of xj in the reference coordinate system if 
ITA = 0. If ITA is greater than zero, then <j>, 9, and ip are 


input. 


1-10 

m 

Direction cosine from x R or <p (deg) 

(Angles read in this 

11-20 

H(Y) 

Direction cosine from y R or 9 (deg) 

order regardless of the 

21-30 

H(2) 

Direction cosine from z R or ip (deg) 

order of rotation(s), ITA) 



Direction cosines of yj in the reference coordinate system if 
ITA = 0. Not used if ITA is greater than zero. 

31-40 

JI(X) 

Direction cosine from x R 


41-50 

JI(Y) 

Direction cosine from y R 


51-60 

JI(Z) 

Direction cosine from z R 


61-65 

IFL 

Must match ITA, Record 1 5, columns 33-35 
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RECORD 17: ANGULAR INTEGRATION CONTROL 

These records enable the program to limit the calculations to the area of interest on 
a conic, a general body of revolution, and a circular plate by limiting the angular surface 
integration. For all other bodies, blank cards should be used. There are as many of 
these cards as there are subshapes and read in order immediately after the subshape 
direction cosines or axes rotation angles (Record 16). 

The choice of PHIM1 and PHIM2 values provides a method of centering the calcu- 
lation nodes such that output data are printed on desired angular locations. The PHIM1 
and PHIM2 angular values are measured in the subshape coordinate system. This does 
not limit the size of the subshape which is always 360 deg. It only limits the areas of 
calculation, thus reducing computer execution time. 

Format: (2E10.6) 


Column I I Variable 

I- 10 PHIM1 

II- 20 PH I M2 


Description 


Angle in degrees measured in the y-z plane from the +y 
axis at which calculations are to be started. 


Angle measured same as PHIM1 at which calculations are 
to be terminated. 
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RECORD 18: FLOWFIELD LIMITS 

This record contains the necessary information to limit the calculations to those areas 
of interest in the flowfield. This scheme is employed in order to make the program 
efficient for the many problem orientations which are possible (see sketch below). Units 
are consistent with IOPT(14), Record 5, column 70. 

Format: (6E10.6) 



I- 10 CUTDAT(I) Radial coordinate defining upper cutoff. 

II- 20 CUTDAT(2) Axial coordinate defining upper cutoff 

21-30 CUTDAT(3) Angle cutoff line makes with horizontal 


31-40 CUTDAT(4) Radial coordinate defining downstream cutoff 
41-50 CUTDAT(5) Axial coordinate defining downstream cutoff 

51-60 CUTDAT(6) Angle cutoff line makes with horizontal 


Ri CUTDAT(1) 

X l CUTDAT(2) 

9i CUTDAT(3) 

R 2 CUTDAT(4) 

X 2 CUTDAT(5) 
d 2 CUTDAT(6) 

0 Normal termination by 
Ri x 2 downstream cutoffs. 

O Normal termination by 
plume boundary. 
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RECORD 19: EXTERNAL FLOWFIELD DIRECTION COSINES OR AXES ROTATION 
ANGLES 

This record contains either the direction cosines of the unit vectors or the three 
axes angles of rotation written in the reference coordinate system for the plume external 
approach flow (EF). (This record is needed only if Record 5, column 38 = 0.) 

Format: (6E1 0.6,15) 

I Column I Variable I Description | 

Direction cosines of x R f in the reference system if ITA = 0, 
columns 63-65. If ITA is greater than zero, then <f>, 9, and ip are 
input 


1-10 

ll(X) 

Direction cosine from x R or <p deg 

(Angles read in this order re- 

11-20 

H(Y) 

Direction cosine from y R or 9 deg 

gardless of the order of rota- 

21-30 

ll(Z) 

Direction cosine from z R or ip deg 

tions), ITA) 



Direction cosines of ysF in the reference coordinate system if 
ITA= 0. Not used if ITA is greater than zero. 

31-40 

JI(X) 

Direction cosine from x R 


41-50 

JI(Y) 

Direction cosine from y R 


51-60 

JI(Z) 

Direction cosine from z R 


63-65 

IFL 

Must match ITA on Record 15, columns 33-35 
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If additional cases are to follow in which a limited number of records are changed, 
Record 5, column 40 = 1, the data are read into the program via a NAMELIST entitled 
CHNGE. The variables contained in the NAMELIST are described in the input cards 
They are as follows: ICONIC, IOPT, TEC, REC, HEADER, ICON, IRUN, PCRAT, TAME, 
IOUT, PAMB, TAMB, ENAMB, GAMAMB, RAMB, CPINF, UNITIC, UNITJC, UNITKC, 
ITU, ITV, ITW, IPLSYM, ISELT, RCI, PC. The program will continue to read cages until 
the last digit in IOPT(8) = 0, Record 5, column 40 or Record 20 ‘END’ is encountered. 


N0TE - Th®variable, ICONIC, contains the information read in on Record 14b, Columns 


NOTE: UNITIC, UNITJC, and UNITKC are input on Record 7. If the angles of rotation 
are being input, all three must be reset through the CHNGE input. 

If IRUN(4) > 0, Record 1 , column 10, Records 2-19 are not used and the impingement 
run previously made and now being input through IRUN(4) data file is updated via 
a NAMELIST entitled INPT2. This allows plotting of the impingement data using the 
following variables in the name list: IOPT, IPLSYM, PRATIO, RICH, PCRAT, ITPLOT 
TUNIT, IRUN, ICON. 


RECORD 20: PROGRAM END 

End of input record. The program will analyze as many cases as the user desires 
for each computer run. The program execution is terminated by end of case card. 

Format: (A4) 


Column | Variable 
1-4 = END 


' ' : ^ ; ^ : ; Description \ y ' : ; . : ^ ' - .■ - ; ; ■ j 

Causes program execution to terminate 
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2.2 Input and Output Files 


The modified PUMP code operates on VAX/VMS and SUN Workstations (UNIX 
operating systems). Before any input or output can take place between a file and a 
program, the file must be connected to a unit number. This task can be accomplished with 
l A h ® ^ of ASSIGN statements in COMMAND files on the VAX system. In keeping with 
ANSI Standard FORTRAN 77 and due to the fact that ASSIGN statements are not easily 
implemented on UNIX, units/files should be connected before use by means of OPEN 
statements. These OPEN statements are located in the FORTRAN file - OPEN. FOR (on 

tape, Ref. [8]). The file names are read using a namelist entitled FILES. The names 
may be input in the data file FILES.DAT. 


The data file, FILES.DAT, is presented in Table 2 along with a description of the input 
variables, and the program magnetic tape assignments are given in Table 3. 

Table 2: Input File to Establish NAMELIST Files 

&FILES 


ORD = 'N' , 

IMP = 'N', 
0RD_IMP = 'N' , 

SORFL = ' Y' , 


ORDER DATA ONLY — INPUT Y 

PERFORM IMPINGEMENT ANALYSIS ONLY — INPUT Y 
ORDER DATA AND PERFORM AN IMPINGEMENT ANALYSIS 
— INPUT Y 

PERFORM AN IMPINGEMENT ANALYSIS USING AN SFPGEN 
PLUME — INPUT Y 


ORD_2 = ' file name for unit 2', IF ORD 

ORD_5 = ' file name for unit 5' , 

ORD_6 “ ' file name for unit 6', 

ORD_8 - ' file name for unit 8', 

ORD _10 = • file name for unit 10', 


' Y # ENTER FILE NAMES FOR 
LOGICAL UNITS 
2,5, 6,8,10 


IMP_2 = ' file name for 
IMP_5 - ' file name for 
IMP_6 = ' file name for 
IMP_8 = ' file name for 


unit 2', IF IMP ' Y' 
unit 5', 
unit 6', 
unit 8', 


ENTER FILE NAMES FOR 
LOGICAL UNITS 
2,5, 6,8 


ord _IMP_ 2 = ' file name for unit 2', IF ORD_IMP = ' Y' 

ORD_IMP _5 = ' file name for unit 5', 

ORD_IMP_6 = ' file name for unit 6', 

ORD_IMP_8 = ' file name for unit 8 ', 

ORD_IMP_l 0 = ' file name for unit 10', 


ENTER FILE NAMES 
FOR LOGICAL UNITS 
2,5,6,8,10 


33 


rsxi t e: o f-n 


RTR 214-01 


Table 2: (Continued) Input File to Establish NAMELIST Files 

S0R_5 = ' file name for unit 5', IF SOR = ' Y' ENTER FILE NAMES FOR 

S0R_6 = ' file name for unit 6' , LOGICAL UNITS 5,6,16 

S0R__16 = 9 file name for unit 16' , 

SEND 

If the user wishes, this option may be removed by commenting out the call open 
statement in main . for and using assign statements to name the different logical units 
on the VAX syste m . 


Table 3: Program Magnetic Tape Assignments 


PUMP 

| Logical Unit 
11 Number 111 

File Function 

IGO = 2 

2 

Limits of ordered flowfield (output) 

(Arrange MOC or 

3 

Scratch file 

RAMP2 data only) 

5 

PLIMP input file (user supplied) 


6 

PLIMP output file 


8 

Ordered flowfield data file (output) 


10 

Flowfield file generated by MOC or 
RAMP2 (input) 

IGO = 0 

2 

Limits of ordered flowfield (input) 

(Perform an impingement 

3 

Scratch file 

analysis ... flowfield has 

5 

PLIMP input file (user supplied) 

been previously 

6 

PLIMP output file 

ordered) 

8 

Ordered flowfield data (input) 


9 

NCAR graphics data file (output) 


11 

Scratch file (used by plot package) 


13 

Coefficients data used in heating analysis 


15 

Scratch file 


18 

Scratch file 


IUNIT 

Variable file unit number on which 
impingement data are to be stored 

IGO = 1 

2 

Limits of ordered flowfield (input) 

(Order flowfield and 

3 

Scratch file 

perform an impingement 

5 

PLIMP input file (user supplied) 

analysis in a single 

6 

PLIMP output file 

execution) 

8 

Ordered flowfield data file (input) 


34 






REIMTECH 


RTR 214-01 


Table 3: (Continued) Program Magnetic Tape Assignments 


PUMP 

Logical Unit 
Number 

File Function 


9 

NCAR graphics data file (output) 


10 

Flowfield file generated by MOC or 
RAMP2 (input) 


11 

Scratch file (used by plot package) 


13 

Coefficients data used in heating analysis 


15 

Scratch file 


18 

Scratch file 


IUNIT 

Variable file unit number on which 
impingement data are to be stored 

IGO = 0; ISORCG = 1 

3 

Scratch file 

or ISORCP = 1 

5 

PUMP input file (user supplied) 

(Perform an impingement 

6 

PUMP output file 

analysis using an 

9 

NCAR graphics data file (output) 

SFPGEN generated 

11 

Scratch file (used by plot package) 

flowfield) 

13 

Coefficients data used in heating analysis 


15 

Scratch file 


16 

SFPGEN plume flowfield file (input) 


18 

Scratch file 


IUNIT 

Variable file unit number on which 
impignement data are to be stored 

IRUN4 > 0 

IUNIT 

Variable file unit number on which 
impingement data are stored 


11 

Scratch file (used by plot package) 


2.3 Pre- and Postprocessing Graphics 

Table 4 shows two examples of run sequences for the PLIMP pre- and postprocessor 
NCAR graphics routine. The first example is the interactive prompts/responses for a 
postprocessing run assuming Record 1 , column 70, NPLOT = 1 and Unit 9 has been 
opened. The second example is a preprocessing run where shading was requested 
and engine coordinates relative to the reference system were given. In both cases the 
user's response is shown boxed. The icon labeled ‘idt’ appears on the screen after the 
program is finished and the picture is ready for display. The mouse should be clicked 
in the ‘display’ box of the icon for viewing the picture. Questions regarding the meaning 
of the prompts are answered in Section 3. 3. 5.3 of Volume I. That section also contains 
a sample input for the preprocessing file using the stand-alone NCAR graphics program 
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without using PUMP. Creation of that file is not interactive and is explained in Section 
3.3.5.2 of Volume I. 


Table 4: NCAR Graphics Run Sequences 


Postprocessing 


Preprocessing 




h ■ 1 1 too L^Z us r/loc* 1/ h i n/b as h, 


dayton> run_n<tr pi lap.? f aortic* 

you WANT TO PLOT AN EXISTING .cga FILE ????(yts/no) 

00 YOU HEED TO CHANGE THE NAFE OP THE PLIN> OUTPUT FUE 

pal 

DOES THE FIIC plcturt.dtt EXIST ? (ytj/no) 
IF YOU ARC HOT SURE PLAY TT SAFE 4ND ENTER yts 
7fti 

oitap.gr tobies 

OO you HEED TO RECOMPILE NCAR FORTRAN yis/tvj 

C53 

GRAPHICS PRCPROCESSIHC > 1 

PIOT PIIK> OUTPUT -> 1 

DIRECTLY PLOT plcturt.dtt > ANY OTHER NO. 


CREATE 1 PLOT plcturt.dtt FROM PtIMP OUTPUT — > 1 
GRAPHICS PREPROCESSOR — > ANY OTHER NO. 


ENTER THE VARIABLE NUWER FOR NOTCH YOU HART CONTOURS 

PRESSURE — > (t) 

TURBULENT HEAT RATE — > C2) 

LAMINAR HEAT RATE — > (3) 

TRANSITIONAL HEAT RATE — > (4) 

FREE MOLECULAR HEAT RATE — > (5) 


fiilkoo 


a ay ton > run_ncar pi 1ao_gr tomes 
DO YOU WANT TO PLOT AN EXISTING .cga FILE ????Cyts/rw) 

no I 

WYOU NEED TO CHANCE THE NAFC OF THE PLIFP OUTPUT FILE 
tnoi 

DOES THE FILE o1cturt.dat EXIST ? (y*f/no) 

IFYOU ARC NOT SURE PLAY IT SAFE AMD ENTER yti 
vew 

pHap.grtohlcs 

DO YOU NEED TO RECCW>ILE NCAR FORTRAN yts/no 

GRAPHICS PREPROCESSING > 1 

PLOT PIIH> OUTPUT -> 1 

OIRECTIY PICT plcturt.dtt — — > ANY OTHER NO. 


REATE 4 PLOT plcturt.dtt FROM PLIW OUTPUT — > 1 
GRAPHICS PREPROCESSOR — > ARY OTHER NO. 


ENTER 39 IF YOU HART THIS SOFTWARE TO SWAOC. 
NOTE THAT SHADING IS VERY COPUTE INTENSIVE AMO 
—JLOWS THE RUN-TIFC CONSIDERABLY. 


COORDINATES RELATIVE TO REFERENCE COORDINATE SYS CX Y 2) 


Display Icon 


„ Copyright (C) 1990 - All Righto Raaarved 
Uhieersity Corporations for Ataospherie Research 
HCAR Tier - U1TCX Version 3.01 


select file 


ORIGINAL PAGE i& 
OF POOR QUALITY 
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Section 3 

SAMPLE PROBLEMS 

This section describes two sample problems and carefully discusses the inputs 
necessary for PUMP to solve them. 

3.1 Sample Case 1 

This section contains a sample case involving the STS-7 mission of June 1 983. It will 
aid the user in familiarization and preparation of input data for the modified PUMP code. 
An RCS plume impingement test was performed to provide force and moment data for 
the Shuttle Pallet Satellite (SPAS-01) due to a single orbiter RCS firing [9]. The SPAS- 
01 was positioned at different orientations relative to the orbiter with the forward upward 
facing RCS engine impinging upon its surface. This test case was previously used as 
an impingement validation case [10]. However a comparison with the experimental data 
will not be performed at this time. The RCS plume flowfield was calculated using the 
RAMP2 computer code. This plume is one of the flowfields generated for the plume 
database described in Volume I. 

The nozzle solution was performed assuming equilibrium/frozen thermochemistry, 
variable oxidizer to fuel ratio distribution. BLIMPJ execution was also required since the 
effects of the boundary layer were important. 

The RAMP2 code cannot perform a plume calculation considering both variable O/F 
ratio distribution and total enthalpy gradients (due to boundary layer effects); therefore, 
the RCS plume flowfield was generated utilizing the ideal gas variable total enthalpy 
option of the code. An SFPGEN plume flowfield was then generated using the RAMP2 
plume file. All the relevant information to perform these tasks is presented in Volume I. 

During the test, the SPAS-01 satellite was positioned at eight different locations 
relative to the orbiter as shown in Fig. 2. The RCS engine was fired while the satellite 
was at each position. There are many options in setting up the coordinate systems in this 
sample problem; Fig. 3 shows the engine and composite (reference) coordinate systems 
chosen. The SPAS-01 structure was divided into 10 subshapes and is presented in 
Fig. 4. The reference coordinate system in Fig. 4 corresponds to the axis and is located 
at the Center of Gravity (CG) of the vehicle. Table 5 lists the PUMP input for this sample 
case and should be referred to while reading the following details for each record. 

Record 1 — Record 1 of the input data concerns the program control information. 
The variable IGO determines the type of calculation required. RAMP2 flowfields have 
to be arranged into a specific format acceptable to PUMP. If the flowfield has been 
previously ordered and an impingement analysis is required, a value of 0 is entered 
for IGO, if not, the data may be arranged into this format and an impingement anal- 
ysis performed in the same execution (IGO - 1) or during a later execution (IGO = 
2). SFPGEN flowfield files are ordered flowfields, therefore, a value of 0 was entered 
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Figure 2: Relative Location of SPAS and Orbiter for Plume Impingement Test 



Figure 3: Orientation of Engine and Composite Coordinate System 
for SPAS/RCS Impingement Calculation 
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Table 5: Input File for Sample Case 1 


0 0 0 0 1 2 000 1 o 

SPAS -01 RCS Plume Impingement Test 


1 

0 3 

10 10 

712.65 

- 199.98 

0.0 

0.0 

180.0 

- 2.00 

0 . 

0 . 


0.1 

- 06250 . 

0.0 

0.0 

0.0 

0.0 

500 . 

500 . 

500 . 

500 . 

500 . 


166 . 

46 . 

166 . 

46 . 

52 . 

46 . 

11.9 

56.53 

56.53 

166 . 

23 . 


2003003000 46 . 

166 . 

00 



2015015002 166 . 

28.0 

00 



2015015002 166 . 

25.6 

00 



2003003002 46 . 

52 . 

00 



2015015002 46 . 

52 . 

00 



4010010002 56.53 

32.5 

3 8 



4010010002 11.9 

56.53 

3 8 9 



4010010002 56.53 

11.9 

3 7 9 



4010010002 32.5 

56.53 

3 6 7 

8 


2003003002 23 . 

166 . 

0 



3.68 

-.18 

14.57 

26.18 

-.18 

.57 

21.67 

-.18 

- 25.43 

3.18 

82.82 

- 11.43 

3.18 

- 83.18 

- 11.43 

11.0 

10.02 

- 37.43 

11.0 

10.02 

- 37.43 

11.0 

- 10.38 

- 37.43 

11.0 

- 10.38 

- 37.43 

5.68 

-.18 

- 37.43 

90.0 

0.0 

90 . 

0.0 

0.0 

0.0 

0.0 

20.6 

0.0 

0.0 

0.0 

90 . 


1 0 000 2 000000 


1 . 

1 . 

1 . 

1.4 

1716 . 

2 . 

2500 . 

1.0 

12 . 

500 . 

500 . 

500 . 

28 . 

166 . 

25.6 

52 . 

56.3 

32.5 

11.9 

32.5 

56.53 


132 

123 

213 

312 

312 

123 

132 

123 

132 

213 


0 0 10 

10 10 
1 231 


1 . 0 . 

2 DFP 

500 . 500 . 


1 
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Table 5: (Continued) Input File for Sample Case 1 


0.0 

0.0 

-90. 

-90. 

-31.2 

23.45 

180.0 

-23.45 

-31.2 

-90. 

25.0 

23.45 

180.0 

-23.45 

25.0 

-90. 

90. 

0.0 

10 blank 

lines (1 

per subshape) for Record 17 

100. 

-100. 

0.0 0.0 100. 


END 


in column 5. The second input variable, IRUN(4) was set to 0 since impingement 
calculations are desired. Therefore IRUN(7) is left blank. Since PUMP can now perform 
two-phase impingement analyses in a single execution LSOLID is no longer required. 
LRAMP is set to 1 as the plume flowfield was generated using the SFPGEN computer 
code. The RAMP2 plume flowfield was generated using the ideal gas with variable total 
enthalpy option, therefore IAF(8) is equal to 2. ISPECI, IDSHAD, and IMTRAN are set 
to 0 respectively as extra output is not required. The next two variables, ISORCG and 
ISORCP, are 1 and 0, respectively (gaseous SFPGEN flowfield.) English units are used 
to prepare the input file (METRIC = 0) and 1 32 column format is desired (FORM = 0). 
The variable NPLOT determines if data for the 3-D NCAR pre/post processing graphics 
package will be output on unit 9. NPLOT was assigned a value of 1 to enable this 
option. If the user wishes to check the input geometry without performing an impingement 
analysis i.e., use the NCAR graphics preprocessor option, ICENT is input as 1 . If the post 
processor option is desired (as in this case), ICENT is assigned a value of 0. ICENT is 
not required if NPLOT is less than 1 . 

Records 2 and 3 — Records 2 and 3 are not required as IGO = 0. 

Note: Sample case 2 in this section demonstrates how to order a RAMP2 flowfield. 

Record 4 — The problem title is input on Record 4 and may be any 72-character 
alphanumeric name. 

Record 5 — This record contains the input options necessary to control program 
execution. The flowfield type is specified by IOPT(1 ) and for this case is input as 1 for 
axisymmetric flow. Continuum and free molecular heating rate calculations are desired, 
thus IOPT(2) is 0. The limit of uniform flow assumptions is given by IOPT(3). It is 
assumed that the elemental areas are in uniform flow (IOPT(3) = 3). The SPAS-01 
satellite was divided into 10 subshapes (IOPT(4) = 10). 1VAPN and IPROP are both 
set to 0 since ISCTAM is 1. IOPT(5) indicates the type of output desired. In this case 
IOPT(5) = 0 as data for each elemental area and full distribution information is required 
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on each subshape. A value of 0 for IOPT(7) indicates the number of species to be input 
for boundary layer calculations. 

Note: The number of species and their mole fractions can be obtained from the 
RAMP2 output file. Not required if using an SFPGEN plume flowfield. 

IOPT(8) is required to scale the desired chamber pressure to that actually used to 
generate the flowfield. No scaling is needed, therefore IOPT(8) is 0. Columns 38 and 39 
determine the method used to calculate the pressure force on the impinged body. The 
plume loading only option is desired. Zero was input in column 38 (therefore, column 
39 is left blank). A zero is also input in column 40 as one case is being executed. The 
type of local impact pressure is specified using IOPT(9). The modified Newtonian option 
is used (IOPT(9) = 2) in this sample case. 

Note: Record 14 provides an opportunity to override this selection for each 
individual subshape. 

A detailed printout of the intermediate steps used in the calculation procedure at 
each subshape may be obtained by inserting into IOPT(10) the number of subshapes for 
which the printout is desired. Since this produces a large amount of data, IOPT(10) was 
input as 0. IOPT(11) is the plot control flag. Columns 51-55 determine the number of 
parameters and the way they would be presented on the plots. No plots are required for 
this case, therefore columns 51-55 are left blank. The number of the computer execution 
is inserted in IOPT(12) and is 1. The units for PUMP are defined by IOPT(13). Inches 
were used in setting up the input file for this sample problem (IOPT(13) = 0). The units 
used to generate the flowfield are feet (IOPT(1 4) = 1). 

IOPT(15) provides the opportunity to print out local intermediate properties during 
the pressure and heating rate calculations. This is set to 0 as printout is not desired. 
Whether the shading option is enabled or not depends on the value of IOPT(16). Shading 
is required so a value of 1 is inserted in column 80. 

Record 6 This record defines the relative positions of the engine and reference 
coordinate systems. The engine position was defined by a vector from the reference 
coordinate system to the engine coordinate system. Thus, column 65 (ITU) is set to 1 . 
Table 6 outlines the relative location of the engine coordinate system with respect to the 
composite coordinate system (SPAS-01 CG) for the eight locations presented in Fig. 2. 
Position 1 was chosen for this case. The X re f axis points towards the engine and results 
in a value of 712.65 for REC(X). The Y Tef axis points away from the engine coordinate 
system, REC(Y) is —199.98. The value for REC(Z) is 0.0. Since no transfer of moments 
is required, RECT(X), RECT(Y), and RECT(Z) are left blank. As a result of this, PUMP 
will use the integrated impingement forces to calculate the resultant moments about the 
reference coordinate system. The engine-reference transformation control flag is input 
in columns 68-70 (ITV). This variable will determine which procedure will be used to 
transform the reference coordinate system to the engine coordinate system. 

Note: The transformation is always from the reference system to the engine system 
regardless of the value of ITU (column 65). 

If the three axes angles of rotation are to be specified on Record 7, the value inserted 
for ITV is the order of this rotation. In this sample problem, ITV = 231, the reference 
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Table 6: Relative Location of Engine Coordinate System with Respect to 
the Composite Coordinate System (SPAS-CG) 


Position Number 

E 

X (In) 

Y (In) 

Z (In) 

1 

712.65 

-199.98 

0.0 

2 

952.65 

-199.98 

0.0 

3 

952.65 

-919.98 

0.0 

4 

1312.65 

-919.98 

0.0 

5 

1312.65 

-199.98 

0.0 

6 

1312.65 

-199.98 

-420.0 

7 

1912.65 

-199.98 

0.0 

8 

2512.65 

-199.98 

0.0 


coordinate is first rotated about the Y axis, then about the Z axis and finally about the 
X axis. 

Record 7 The numerical values of the three axes angles of rotation are inserted 
on this record. The X axis rotation is 0.0 deg, the Y and Z axes rotations are 180 deg 
and -2.0 deg, respectively. Figure 5 demonstrates these rotations in pictorial form. The 
values are input in order of X, Y, Z regardless of the rotational order specified by ITV 
Record 6. 

Record 8 — - This record is used to identify the parameters for which plots are desired. 
It is input only if ITPLOT > 0 (IOPT(1 1 )) on Record 5. In this case, this record is not input. 

Record 9 — Record 9 is the species concentration record. It specifies the mole 
fractions of the flowfield gas composition when the species are not in the input flowfield 
file (Unit 8). It is input when IOPT(7) > 0 and IOPT(2) < 2 on Record 5. This record was 
not used in this case because the input flowfield was generated by SFPGEN. 

Record 10 Engine references and parameters for non-dimensionalizing are given 
on Record 10. XTE is the transfer distance from the engine gimbal to the exit plane. This 
is primarily used when the orientation of the engine coordinate system with respect to 
the reference system is about the engine gimbal point. The units must be consistent with 
the units of the ordered flowfield (IOPT(14)). Since no transfer was desired, XTE was 
input as 0.0. XMIN is the minimum distance from the engine reference system at which 
impingement calculations are to start. This variable is used to reduce computer execution 
time. If it is known that no impingement will take place for some distance from the nozzle 
exit, this distance may be inserted thus eliminating unnecessary calculations. A value 
of 0.0 was used for XMIN. The engine scale factor (REFERN) can be used to convert 
a non-dimensionalized plume to dimensional units or to scale the plume. A value of 1 .0 
was input as neither option was desired. The diameter of the nozzle at the exit plane, 
RICH, can be used to non-dimensionalize the program output. A value of 1 .0 was input 
to avoid non-dimensionalization. The pressure output can be non-dimensionaiized by 
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Engine Coordinate System 



Rotation of 180° about y& Axw Rotation of - 2° about z Axis 


Figure 5: Rotation of the Reference Coordinate System into Alignment 
with the Engine Coordinate System 

the variable PC (engine chamber pressure). PC was assigned a value of 1 .0 since non- 
dimensionalization was not desired. PCRAT is used to change the heating rate values, 
forces and moments to correspond to a desired plume chamber pressure instead of the 
actual chamber pressure used to generate the plume flowfield. PCRAT was input as 1 .0 
as no change is needed. DELPHI allows impingement calculations to be started away 
from the stagnation line; 0.0 deg was input. 

Record 11 — Conditions external to the flowfield are input on this record. The 
ambient static pressure, PAMB, is 0.1E-06 psf. The ambient static temperature, TAMB, 
is 250.0 R. These conditions correspond to an altitude of approximately 300,000 ft. 
EMAMB is 0.0. The ratio of specific heats (GAMAMB) and the gas constant (RAMB) 
for the ambient air are 1.4 and 1716.0 ft^/sec^/R. The Newtonian pressure coefficient, 
CPN, was input as 2.0. 

Record 12 — Record 12 contains the heating rate reference parameters. A value of 
0.0 BTU/lbm was input for CPS (the specific heat for solid particles) since there are no 
solid particles in this plume. The trajectory time, TAME, is input if surface temperature 
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calculations are to be made. TAME was set to 0.0 seconds. The axial distance (XSTAG) 
from the subshape origin to the stagnation point is 0.0 in. The temperature (TD) of the 
gas in the boundary layer at the stagnation point at onset of dissociation is assumed to 
be 2500 R. The reference radius, RS was input as 1.0 in.; therefore, FUDGE was input 
as 12.0. The nose shape (SHAP) of the vehicle is a rectangular flat plate (2DFP). 

Records 13a, 13b, and 13c — Record 1 3a, the subshape surface temperatures, are 
input only if heating rate and/or impingement gaseous contaminant calculations are to be 
performed. The wall temperature, TWAL, for each subshape was set to 500 R. Records 
13b and 13c are required only if the user has selected the deposition model (ISCTAM 
= 2) and has indicated that there is additional surface evaporation rate data for species 
that are not already included in the internal program data. These records were not used. 

Records 14a, 14b, 14c, and 14d — Record 14a contains the body reference 
dimensions, REFL and REFD (reference length and diameter), used to determine the 
local flow regime for the N subshapes. The properties for three subshapes may be listed 
on each line. In this sample problem, ten subshapes were listed. Records 14b, 14c and 
14d are input in succession for each subshape. 

The type of subshape (conic, rectangular plate, circular plate, etc.) are outlined on 
Record 14b. The SPAS-01 structure was defined using 10 subshapes. Six of these 
subshapes are rectangular flat plates, the remaining four are triangular subshapes The 
10 subshapes are outlined in Fig. 4 of the SPAS-01 geometry. The first subshape was 
a flat plate of dimensions 46.0 in. and 166.0 in. Therefore, IDSURF in column 1 is 
2. The surface integration controls were selected with three divisions in the Z direction 
(NPHI = 3) and three divisions in the Y direction (NX = 3). The values chosen for the 
integration controls are dependent on the problem. The flat plate is composed of 9 (3 x 3) 
elemental areas as a result of assigning these values to NPHI and NX. The total number 
of elemental areas per subshape will effect computation time, therefore unnecessary 
calculations should be avoided. 

The type of impingement pressure calculation (if different from the option specified 
by IOPT(9) on Record 5) is recorded in column 8. For this case column 8 was input 
as 0, thus the option previously chosen will be used for this subshape i.e., modified 
Newtonian pressure distribution. Column 9 indicates the heating rate option for each 
subshape. Since heating rate calculations are not required for this subshape, columns 
9 and 10 are both input as 0. The reasoning behind this decision is straightforward. A 
visual inspection of Fig. 4 shows that this subshape is shaded from the flow by all other 
subshapes when the engine is located in position 1. The coefficients of the subshape 
equations are input in columns 11 - 60. The coefficients are defined in the subshape 
coordinate systems. The value of A is 46.0 in. and B is 166.0 in. Since this subshape 
is a flat plate; there are no values for C, D, E, XN or XB. 

The shading option [IOPT(16)] was flagged on Record 5, therefore Record 14c is 
required. Zero is input in column 2 of this record. This informs PUMP that no other 
subshape will be checked to see if it is shading this one. Although each elemental area 
comprising this subshape will be checked to determine whether or not it is pointing in 
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the direction of the flow. If they are pointing away from the flow (as in this case) the 
subshape shades itself. 

The final record (in this set) for subshape 1 is a surface temperature identification 
record. Record 14d was not input for this sample problem. 

Subshape 2 is also a flat plate of dimensions 1 66.0 in. and 28.0 in. The surface 
integration controls were selected with 15 divisions in the Z direction (NPHI = 15) and 15 
divisions in the Y direction (NX = 15). Column 9 was assigned a value of 0 to indicate 
that heating rate calculations (using the Eckert Reference Enthalpy method, IHT = 2) 
are to be performed on all elemental areas of this subshape. Again only coefficients 
A and B are input (A = 166.0, B = 28.0). Subshape 2 is facing the flow direction and 
the flow impinging upon it is not impeded by any other subshape; therefore, column 2 
of Record 14c is 0. Record 14d was not used for this subshape. Records 14b, 14c 
and 14d are similar for subshapes 3, 4 and 5 as they are also flat plates and are not 
shaded by other subshapes. 

Subshape 6 is a triangular plate (IDSURF = 4). NPHI and NX were both assigned 
identical values of 10. Therefore the plate was divided into 55 elemental areas. The 
method used to calculate the impingement pressure remains as input on Record 5. 
Eckert Reference Enthalpy was chosen as the desired heat transfer method. A and B 
for subshape 6 were entered as 56.53 in. and 32.5 in., respectively. There are no values 
for C, D, E, XN or XB. 

The flow to subshape 6 may be blocked by subshapes 3, and 8; therefore, 3 was 
inserted in column 6 of Record 1 4c, and 8 in column 4. This specifies that subshapes 
3 and 8 will be checked to see if they block flow impinging on subshape 6. Record 14d 
was not used for this subshape. 

Subshapes 7, 8 and 9 have similar inputs regarding record 14b. Impingement on 
subshape 7 may be blocked by subshapes 3, 8 and 9. Therefore, Record 14c for 
subshape 7 is input as 3 in column 2, 8 in column 4, and 9 in column 6. Similar 
observation shows that subshape 8 may be impeded by subshapes 3, 7 and 9, while 
subshape 8 is screened by 3, 6, 7 and 8. 

Note: If the vehicle were comprised of subshapes generated using the conic, 
circular plate, body of revolution, or trapezoid equations, the coefficients 
of these subshapes would be input according to the instructions outlined 
in this input guide (Record 14b). 

Record 15 — There are 10 subshape position vector records (Record 15) input as 
a group immediately after Record 14. The origin of each subshape is located by a 
vector (RCI) from the reference coordinate system to the subshape origin. The units are 
consistent with IOPT(13) and these records are input in the same order as the subshape 
geometry data (Record 14). For subshape 1 of the sample problem, the X component 
of the position vector was input as 3.68. The Y component was input as -0.18 and 
the Z component as 14.57. All of these distances are in inches. The position vector 
components indicate the location of the subshape origin. A system of axis rotations 
must be chosen to define the subshape axis positions and inclinations. A reference- 
subshape system flag (ITA) is located in columns 33-35 of this record. If column 33 
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is 0, the direction cosines defining the orientation of the two coordinate systems are 
input. If ITA is greater than 0, the input is the order of rotation needed to align the 
axes of the reference system to that of the subshape. It is important at this moment 
to discuss the significance of subshape coordinate placement. Continuum heating rate 
results are highly dependent on the proper placement of the subshape origin. Referring 
to the different subshape descriptions on pages 24-25: 

1 . The X coordinate must be placed along the body centerline in the general direction 
of the flow for conic bodies and bodies of revolution. 

2. For rectangular flat plates, the coordinate system is located at the geometrical center. 
It is located at the vertex of the right angle for right triangular plates. For both, the X 
coordinate must be placed perpendicular to the surface and must be pointing toward 
the general direction of flow. 

3. Circular plates have the coordinate system at the geometric center. The X axis 
is normal to the plane of the plate. Orientation of the Y, Z axes directions are 
established by the user. 

4. The trapezoid subshape coordinate system must be correctly orientated as shown 
on page 25 of this guide. The X axis is normal to the surface. The Y and Z axis 
directions are defined by the user. 

Record 16 — This record contains the subshape direction cosines or the axes angles 
of rotation. There is one of these records for each subshape. In this problem, the rotation 
angles have been input for all ten subshapes. It is important to note that no matter in 
which order the individual rotations are to occur (specified by ITA on Record 15), the 
angles of rotation are input always in the following order: 

1. Rotation about the X axis (columns 1-10) 

2. Rotation about the Y axis (columns 11-20) 

3. Rotation about the Z axis (columns 21-30) 

The rotations for subshape 1 are as follows: the rotation about the X axis was input as 
90.0 deg; the Y axis rotation was 0.0 deg; the final rotation is 90.0 deg about the Z axis. 
Since ITA was 132, the reference coordinate axes are first rotated 90.0 deg about the X 
axis. The second rotation is also 90.0 deg but about the Z axis this time. The reference 
coordinate system is now identically oriented with the axes of this subshape. Therefore, 
the Y axis rotation was specified as 0.0 deg since no more rotations are necessary. The 
input for subshape 2 was 0.0 in columns 1-10, 11-20, and 21-30. This implies that the 
reference coordinate axes are aligned with those of subshape 2. ITA (column 33-35 on 
Record 15) was input regardless of this fact. The angles of rotation for the remaining 
six subshapes were input according to the order specified on Record 14. The numerical 
values may be identified from Table 5. 

Record 17 — This is the angular integration control record. This enables PUMP 
to limit the calculations to the area of interest on conical subshapes, general bodies of 
revolution and circular plates. PHIM1 and PH I M2 are defined in the subshape coordinate 
system. PUMP calculations will only be performed between these limits. The angle is 
input in degrees measured in the Y-Z plane from the +Y axis. PHIM1 specifies the 
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angular location at which calculations ara to bagin, PHIM2 specifies where calculations 
are terminated. There are as many of these records as there are subshapes. Since the 
SPAS-01 was described using rectangular flat plates and right triangular plates, 10 blank 
records are inserted into the input file. It is important to remember that, regardless of 
the fact that only impingement parameters in the area between PHIM1 and PHIM2 are 
calculated, the entire 360 deg is used for shading calculation. 

Record 18 — The necessary information to limit the calculations to the areas of 
interest in the flowfield is input on Record 1 8. The calculation limits are used to create 
a box either about the entire flowfield (as in this case) or about a particular portion of 
the flowfield. PUMP checks for possible impingement only within the specified limits of 
this box. This option is important, particularly for high altitude plumes (large plumes). 
If a subshape is located far downstream from the nozzle exit, these limits would allow 
PUMP to avoid the search through unnecessary portions of the flowfield for impingement. 
For the sample problem the following input was used: columns 1-10 = 100.0, columns 
11-20 = -100.0, columns 21-30 = 0.0, columns 31-40 = 0.0, columns 41-50 = 100.0, 
and columns 41-60 - 90.0. These limits are similar to those used to calculate the 
RAMP2 plume, the only difference being the downstream cutoff limit (X2) was 25.0 ft. 
for the RAMP2 calculation. Since the SPAS-01 and the RCS engine are 59.8 ft. apart 
at position 1, the RAMP2 flowfield is too short for this analysis. The SFPGEN flowfield 
was used. (Refer to Volume I about SFPGEN flowfield length limitations). 

Record 19 — This record describes the external flow coordinate system, but is not 
used in this case since column 38 of Record 5 is equal to zero. 

Record 20 Record 20 was input in columns 1-3 as END. This is used to identify 
the completion of the input data. 

This completes the input file for this sample problem. Selected sample pages of the 
output are in Appendix A. 


3.2 Sample Case 2 

In order to perform an impingement calculation using a RAMP2 flowfield, the flowfield 
file (unit 3) must be ordered into a format acceptable to PUMP. The PUMP code must 
be executed using the first 3 records in Table 7. The PUMP output files, unit 2 and unit 
8 can then be saved and used as input for all subsequent PLIMP executions involving 
this particular flowfield. 

Table 7: Sample Case 2 Input 


20 0012 000 000000 
121 0 0 0 346 2 0 0 

0.0 .775 1.0 

SPAS-01 RCS Plume Impingement Teat - file to create an ordered flowfield 
___ of RCSVOF.UN3 
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Record 1 — Record 1 of the input data concerns the program control information. 
The variable IGO determines the type of calculation required. It is set to 2. This informs 
PUMP to arrange the RAMP2 flowfield file, but not to perform an impingement analysis. 
The second input variable, IRUN(4), was set to 0, with IRUN(7) left blank. LRAMP is set 
to 1 (RAMP2 flowfield file). The RAMP2 plume flowfield was generated using the ideal 
gas with variable total enthalpy option, therefore IAF(8) is equal to 2. ISPECI, IDSHAD, 
and IMTRAN are set to 0 respectively as extra output is not required. The next two 
variables, ISORCG and ISORCP, are set to 0. English units are used to prepare the 
input file (METRIC = 0) and 132 column format is desired (FORM = 0). The variable 
NPLOT determines if data for the 3-D NCAR pre/post processing graphics package will 
be output on unit 9. This option is not required, NPLOT was input as 0. ICENT is not 
required if NPLOT is less than 1 , therefore ICENT is 0. 

Record 2 — Since IGO was set to 2, Record 3 must be input to control the ordering 
of the plume flowfield. ISTART on column 5 determines the characteristic line number 
in the flowfield at which the data ordering process will begin. This corresponds to 
the characteristic line number of the nozzle exit plane and the value is determined by 
examining the RAMP2 output file (unit 6). ISTART is input as 121. ISIGN concerns 
MOC generated flowfields and is set to 0 for this case. Points may be deleted from 
the flowfield by setting IDEL to a value greater than zero. This option was not chosen. 
No intermediate data is to be printed, so IPRINT is 0. ITERM is the characteristic line 
number at which the ordering process is to terminate. This value is determined from 
the RAMP2 output file. It is input as 346 for this flowfield. The cutoff limits specified in 
the RAMP2 input file (unit 5) are to be used to check if a point is within a prescribed 
boundary; therefore, ISEND is 2. ISKIP is set to zero since ISEND = 2. No deletion of 
characteristic lines is required, ISKPL on column 45 is 0. 

Record 3 — Record 3 is used to make any desired coordinate system translations. 
The variables RREF and XREF were used to move the origin of the flowfield from the 
nozzle throat to the nozzle exit. RREF was input as 0.0 ft and XREF as 0.775 ft. The 
units are consistent with the plume dimensions. 

Since IDEL (Record 2, column 20) is zero, DELETE is left blank. No points will be 
deleted from the flowfield. The last variable on this record, DIAM is set to 1 .0 (no scaling). 

No more data is required to order RAMP2 flowfield files. Although (as in this case) 
additional information, such as the problem title, may be input as a reminder. 

Note: If IGO = 1 (Record 1, column 5) l.e., arrange the data and perform an impinge- 
ment analysis, Records 4-20 would be input (according to the particular 
problem) Immediately following the data in this file. 

Record 20 — END is put in columns 1-3 to end the run. 

This completes the input for this sample case. Selected pages of the output are 
shown in Appendix B. 
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Appendix A 

SAMPLE CASE 1 OUTPUT 

This appendix presents selected pages from the complete output file (UNIT 6) for 
sample case 1. The program output is organized sot hat the initial pages contain the 
input data. Impingement data for each subshape is then listed in input order. The 
calculations were performed in English units. The complete output is 21 2 pages. 
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SHEAR FORCE COMPONENTS TORQUE COMPONENTS 

FORCE (X) FORCE (Y) FORCE (Z) TORQUE (X) TORQUE (Y) TORQUE (Z) 

I.64117E+00 0. 35381E+00 -0.31120E+00 0.10487E+00 0.16130E+01 -0.21421E+01 
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Appendix B 

SAMPLE CASE 2 OUTPUT 

The following section lists the complete output file (UNIT 6) for sample case 2. 
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